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Abstract 

Nonlocal order parameters for deconfincmcnt, such as the entanglement entropy and 
Wilson loops, depend on spatial surfaces E. These observables are given holographically 
by the area of a certain bulk spatial surface T ending on E. At finite charge density 
it is natural to consider the electric flux through the bulk surface T in addition to its 
area. We show that this flux provides a refined order parameter that can distinguish 
'fractionalized' phases, with charged horizons, from what we term 'cohesive' phases, with 
charged matter in the bulk. Fractionalization leads to a volume law for the flux through 
the surface, the flux for deconfined but cohesive phases is between a boundary and a 
volume law, while finite density confined phases have vanishing flux through the surface. 
We suggest two possible field theoretical interpretations for this order parameter. The 
first is as information extracted from the large TV reduced density matrix associated to 
E. The second is as surface operators dual to polarized bulk 'D-branes', carrying an 
electric dipole moment. 



1 Holographic phases of matter 



Certain large N gauge theories admit a hole-graphically dual description in terms of a theory 
of classical gravity in one higher spatial dimension flj. Many of the most interesting questions 
one would like to ask of such theories concern the gauge field dynamics. For instance, 
one would like to gain a dual geometrical understanding of confinement. This was largely 
achieved shortly after the discovery of holography [2j, as we will briefly review. 

If the large N gauge theory has more structure, one can ask more refined questions. In 
particular, if the gauge theory has a global £7(1) symmetry, then one can place the theory 
at a nonzero chemical potential /i, which will typically induce a charge density (J*). In this 
finite density state, one can ask whether the charge carriers are "confined" or "deconfined". 
As we shall elaborate throughout this paper, this is not the same notion of confinement 
and deconfinement as we have used in the previous paragraph, which referred to the charge- 
neutral glue sector. New concepts are required. Phases in which the charge carriers are 
deconfined are known as fractionalized phases of matter Q We will introduce the notion of 
a "cohesive" phase to describe the opposite of a fractionalized phase, i.e. when the charge 
carriers are confined j£] 

Various (nonlocal) field theoretic order parameters exist for deconfinement. Confining 
and deconfined phases can be distinguished using Polyakov loops, Wilson loops, and the 
entanglement entropy. All of these have beautifully simple realizations in holography as 
certain surfaces in the bulk geometry (21 El SI E] • Deconfined phases are characterized 
by an event horizon in the bulk, while confined phases are dual to geometries that are 
effectively cut off at some IR scale. Similarly sharp order parameters do not yet exist for 
fractionalization, as distinct from deconfinement. The objective of this paper is to start to 
remedy this situation from a holographic perspective. 

It has been proposed that fractionalization can be diagnosed by a mismatch in the 
Luttinger count in phases that do not break the U(l) symmetry spontaneously H2 El IE]- 
This suggestion led directly to the observation that finite density holographic duals are 
naturally classified by whether the asymptotic electric flux emanates from behind an event 
horizon or is carried by charged matter in the bulk, e.g. [9]. For fermionic bulk matter, the 
flux from behind the horizon indeed appears as a deficit in the Luttinger count. An elegant 

1 "Fractionalized" is often used synonymously to "deconfined" in a condensed matter context, and also 

used to describe a class of gapped phases, but we will use the term more restrictively. 

2 The term "mesonic" has previously been used by one of the authors to describe the opposite of frac- 

tioanlization, but ultimately this may have misleading associations. 
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general proof of this statement was obtained in [10J, building on previous partial results 
[TTJ Q2J [13]. Furthermore, this picture has lead to the interesting suggestion that in the 
superfluid case, where the U(l) is spontaneously broken, fractionalization could be detected 
as a deficit in the transverse Magnus force acting on a moving superfluid vortex |10] , 

It is desirable, however, to have a clean field-theoretic order parameter characterizing 
fractionalization, in the spirit of the various nonlocal order parameters defining deconfine- 
ment. In a holographic context, following the logic of the previous paragraph, this translates 
into finding a quantity that can be measured at the conformal boundary of the dual space- 
time, and that detects whether the asymptotic electric flux originates from behind a horizon 
or not. In this paper we will construct two, closely related, such quantities. We do not 
know at this point precisely what the objects we consider correspond to in the boundary 
quantum field theory. This is clearly an important question. We will argue that they corre- 
spond in the first case to a refinement of the entanglement entropy that can be defined in 
systems at a finite charge density, and in the second case to operators describing branes in 
the bulk carrying a dipole moment. Let us first recall relevant established results concerning 
confinement in the neutral 'glue' sector. 

1.1 Neutral sector: Confinement and deconfinement 

We will briefly review how confinement and deconfinement appear holographically. We will 
focus on the entanglement entropy. The entanglement entropy is associated with a spatial 
hypersurface S^— l in the d+1 dimensional field theory. 

The holographic prescription for calculating entanglement entropy has been proposed by 
Ryu and Takayanagi [5j: given a boundary hypersurface E, find the minimal surface T in 
the bulk that ends on S, i.e. that satisfies X = dT. To leading order in the bulk semiclassical 
limit, the entanglement entropy is then 

where Af is the area of the bulk minimal surface, and Gat is the bulk Newton's constant. 
The area will be divergent because the minimal surface reaches all the way to the asymptotic 
boundary. It is natural to cut off this divergence at a UV scale, where it provides the short 
distance "boundary law" of the entanglement entropy. 

Holographic computations of the entanglement entropy appear rather similar to those 
of spatial Wilson loops [3j 0]. In 2 + 1 field theory dimensions in particular, both quantities 
depend on a spatial curve Si. However, there can be subtle differences. The entanglement 
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entropy is given by the area of a minimal bulk spatial hypersurface in the Einstein frame 
metric. Areas of hypersurfaces in the Einstein frame have the important property that they 
are preserved under dimensional reduction. This allows the entanglement entropy to be 
computed without depending sensitively on the bulk UV completion. Spatial Wilson loops, 
in contrast, are computed by the bulk on-shell fundamental string worldsheet action in a full 
ten dimensional string theory background. This is consistent perhaps with the fact that the 
Wilson loop refers to a specific gauge-theoretic operator, whereas the entanglement entropy 
can be defined independently of the operator content of the quantum field theory. 

The question of interest is then how scales with the linear dimension R of the bound- 
ary region, in the limit of large R. There is an unfortunate linguistic clash here between 
the terms commonly employed for Wilson loops and entanglement entropy, with "area law" 
taken to mean different things. We will refer to a "boundary law" as a scaling 

A r ~ R*- 1 . (2) 

In contrast, a "volume law" will be 

A T ~R d . (3) 

We have already noted that the entanglement entropy universally has a UV sensitive bound- 
ary law contribution due to short distance correlations across the surface E. 

Geometries dual to confining theories typically have the property that the effective (d+2)- 
dimensional spacetime metric collapses in the far IR. In many circumstances this can be 
understood as an internal dimension capping off [2]. The collapse reduces the number of 
low energy degrees of freedom, typically leading to a mass gap. In terms of the hypersur- 
faces needed to compute the entanglement entropy, the collapse of the geometry allows the 
hypersurfaces to terminate in the far IR, as their cross sectional areas vanish. 

Consider the case where the hypersurface E consists of two parallel infinite spatial hyper- 
planes separated by a distance L. We can think of this as the limit of a rectangle where all 
lengths except one have been taken large. With a confining gravity dual, there are then two 
candidate minimal surfaces T. One connects the two hyperplanes, while the other is discon- 
nected, connecting each hyperplane separately to the far IR region where the bulk surfaces 
end. It has been found in several examples that at large separations L, the disconnected 
surfaces have a lower area and therefore determine the entanglement entropy |144 [T5| I16j . 
This leads to an entanglement entropy at large L with a pure boundary form 

S E {L) ~ + Vol (E) • const. , (4) 
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where e is a UV cutoff. In particular, the entanglement entropy is independent, at leading 
order, of the separation L. This process is illustrated in figure [T] below. 




IR 




L 



IR 




L 




Figure 1: Left, minimal bulk hypersurfaces in a confining geometry. At large boundary 
separations L the surface is disconnected and leads to a boundary law for the entanglement 
entropy. Right, minimal hypersurfaces in a finite temperature black hole geometry. At large 
separations L the surface remains connected and leads to a volume law. In both cases, the 
surfaces for small boundary separations are connected. 

Geometries dual to deconfined phases, in contrast, are typically characterized by an 
event horizon in the IR. Horizons can be scale invariant, such as the AdS Poincare hori- 
zon or the various extremal horizons discussed in detail below, or can be finite temperature 
Schwarzschild-like horizons. Zero temperature horizons are considered in the following para- 
graph. For finite temperature horizons, similarly to the case of confining geometries, there 
is an IR scale at which the geometry ends. However, the geometry close to the horizon is 
distinct: minimal spatial surfaces are unable to cross the horizon [5j[T6l[T7]. Instead, as the 
boundary hyperplanes we considered in the previous paragraph are taken far apart, the bulk 
surfaces droop down close to the horizon and then extend along the horizon. This leads to 
a volume law of the entanglement entropy at large separation 



The second term here is the volume scaling term and is essentially just the thermal entropy 



For extremal horizons, while the spatial bulk surfaces are still unable to cross the horizon, 
the emergence of a scaling symmetry in the far IR typically results in a bulk surface that does 
not lead to a volume law. The IR correction to the UV sensitive boundary scaling is a term 




(5) 



of the system. This behavior is also illustrated in figure [T] 



4 



in the entanglement entropy that vanishes at large separations, scaling like L 9- ^" 1 ), where 
9 < d — 1 is the hyperscaling violation exponent |18[ I19j . We will derive this fact below. If 
9 = d — 1 one finds a logarithmic enhancement to the boundary law |18[ [T9] . One interesting 
exceptional case is the extremal Reissner-Nordstrom horizon, where a volume scaling term 
as in ([5]) does appear due to the ground state entropy density and corresponding finite size 
horizon |20| . Other exceptional cases are locally critical theories without a ground state 
entropy density, recently discussed in |21j . In these cases, at large separation the bulk 
hypersurface appears to be disconnected, leading to an entanglement entropy of the form 
Q, as for confining geometries. 

1.2 Charged sector: Cohesion and fractionalization 

The previous subsection reviewed how the entanglement entropy distinguishes between var- 
ious confining and deconfined holographic phases. The essential point is that the bulk 
minimal surfaces ending on well separated parallel boundary hyperplanes probe the far IR 
spacetime geometry. These have qualitatively different behavior depending on whether the 
IR geometry collapses or ends at a horizon. We would like a refinement of the entangle- 
ment entropy that can determine whether there is flux emanating from behind a horizon, in 
addition to detecting the horizon itself. 

In figure [2] below we show the various classes of bulk backgrounds that we would like 
to differentiate. We are primarily interested in zero temperature solutions. There are at 
least four possibilities. If there is a horizon in the spacetime, then all or some of the electric 
flux may emanate from behind the horizon (we will include certain null singularities in our 
notion of horizons). These are the fractionalized and partially fractionalized cases. If, in 
contrast, all flux originates from charged matter in the spacetime, then there may either be a 
neutral horizon in the far IR, or the geometry may confine. These are deconfined cohesive 
and confined cohesive phases, respectively. One may find in addition spacetimes in which 
the charge emanates from a singularity, but with a geometry that is confining according 
to the usual criteria |22j . Known examples have apparent pathologies such as hyperscaling 
exponent 9 > d—1. For this reason we have not included a picture corresponding to confined 
but fractionalized phases. It is of interest to see if such phases can exist. 

All of the possibilities shown in figure [2] have been constructed. Fully fractionalized 
phases include extremal black holes with and without dilaton fields |23[ |2~H [25j [26l [271 EH] . 
Partially fractionalized backgrounds were constructed in |29[ 130] , Backgrounds with de- 
confined neutral horizons and cohesive charge carriers include electron stars |31[ [32] I33j . 
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Figure 2: From left to right: (i) all flux emanates from horizon, (ii) a fraction of the flux 
emanates from a horizon, a fraction is from charged bulk fields, (iii) all flux from bulk fields 
outside a neutral horizon, (iv) all flux from bulk fields in a confining geometry. 



the ground states of holographic superconductors |34l 135] and the topologically charged 
spacetimes of |30| , Confined cohesive phases include electron stars and holographic super- 
conductors in confining spacetimes |36| 1371 |38j |39] . 

The entanglement entropy alone is unable to distinguish between charged and neutral 
horizons. This is clear because it is possible to have the same near-horizon geometry with 
and without flux emanating from behind the horizon (compare for instance |26| and [32j). 
The area A-p of a bulk spatial hypersurface Y is the simplest invariant one can associate 
to the surface. However, given a bulk theory with a Maxwell field, a second equally simple 
invariant is the conserved electric flux 3>p through the hypersurface. The origin of the electric 
flux in the bulk is precisely the question we are trying to address. Therefore a flux-sensitive 
observable is likely to be of use. This motivates consideration of the following 'deformed' 
holographic entanglement entropy 

^ = j^+7*r. (6) 

The idea is that we will again associate a bulk spatial hypersurface Y to a boundary spatial 
hypersurface £ such that dY = S. The deformed entropy (|6| is to be evaluated on this 
bulk surface Y. In the remainder we will discuss exactly which bulk hypersurface Y should 
be used in the expression for the deformed entropy. We will consider two scenarios. In the 
first, we take the surface Y to be the same Ryu-Takayanagi minimal surface as is used for 
the entanglement entropy. In this case, the flux $p can be thought of as additional data 
associated to the entanglement entropy. In the second scenario, we extremize the deformed 
entropy itself. Here <E>r is to be thought of as an electric dipole moment that has been 
induced on a D-brane worldvolume theory. Such polarized D-branes would presumably be 
dual to a class of surface operators. 
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2 A hypersurface order parameter 



For any bulk spatial hypersurface T, the flux $r is easily evaluated in fully fractionalized 
cases. The absence of explicit charged sources, combined with Gauss's law and assuming 
for the moment that the surface does not fall through the horizon (we recalled above that 
minimal spatial hypersurfaces indeed cannot), leads to the conclusion that 



Thus $r obeys a volume law in these cases. Here p is the charge density of the boundary 
field theory and we have again taken the asymptotic boundary E of T to be two very large 
parallel spatial hyperplanes separated by a distance L. Recall that the field theory charge 
density is given by the electric flux at the asymptotic boundary (9J . The result (JTl) is proven 



Figure 3: By Gauss's law, in the absence of bulk charges, the electric flux through the 
surface T is equal to the flux at the asymptotic boundary in the region bounded by the 
surface X = dT. The total asymptotic flux is pL ■ Vol(S). 

By <3?r we will always refer to the conserved flux. In the presence of bulk dilaton fields, 
this will not be simply the integral of the Maxwell field strength. Gauss's law, using if 
necessary a modified flux due to dilatons |29[ HO] , plays a key role in holographic derivations 
of the Luttinger theorem [10, 13j. It seems natural that it should also feature prominently 
in our discussion. 

Another situation that is easy to describe is when there are disconnected bulk surfaces 
that fall straight into the IR and terminate there, such as the minimal surfaces in confining 




(7) 



in figure [3] below. 
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geometries that we reviewed above. In this case, independently of whether the flux originates 
from a horizon or from charged bulk matter, there is no flux through the bulk hypersurface 
and so 

$r = . (8) 

Here we are making the assumption that the bulk flux is purely radial and that the bulk 
surfaces fall straight into the IR without bending and picking up some flux. 

The cases that remain to be considered are those in which at least part of the flux is 
sourced by charged matter in the bulk and where the hypersurface T is not of the discon- 
nected type described in the previous paragraph. These cases are more sensitive to the 
precise form of the surface involved. In the following two subsections we will consider two 
distinct natural choices for the bulk hypersurface. As we discuss in the final section, these 
correspond to two distinct dual field theoretic interpretations of the flux. 

2.1 Electric flux through minimal hypersurfaces in the bulk 

Minimal spatial hypersurfaces in the bulk determine the entanglement entropy of boundary 
regions according to the Ryu-Takayanagi proposal reviewed above. The simplest choice of 
hypersurface T on which to evaluate the deformed entropy would thus seem to be these min- 
imal hypersurfaces. From this perspective, the deformed entropy is not really being taken 
seriously as a new observable; rather we are associating a new quantity, the flux through 
the bulk hypersurface, to the data that arises in computing the entanglement entropy holo- 
graphically. 

We have already computed the flux through minimal hypersurfaces in fully fractionalized 
([7]) and in confining phases To discuss partially fractionalized and deconfined cohesive 
phases, we need to find the explicit form of the hypersurface in the far IR. This depends 
of course on the spatial IR geometry of the bulk spacetime. We will consider a large class 
of zero temperature scaling geometries that have arisen in recent holographic studies as the 
near horizon geometries of extremal black holes, both with |31[ |32| l29j [30] and without 
|25[ [26[ [27] charged bulk matter. The IR spatial geometry takes the form 

J* 2 = dr2 + (a\ 

"^spatial r 2(d-20)/(d-8) ^ r 2 ' W 

Here the IR limit is understood to mean r — > oo. The boundary spatial dimension is d and 
9 is the hyperscaling violation exponent |19| [28] . The dynamical critical exponent z will not 
appear in our discussion. 



8 



We wish to consider hypersurfaces that are extended along d — 1 of the boundary spatial 
dimensions, forming a large rectangle with area Vol (S), and that form a curve x(r) in the 
remaining dimension. In that case the area of the hypersurface is given by 

A = Vol(S) J *V r 2e/(d-e) + i ,2_ ( 10 ) 

A minimal hypersurface T is then easily found to solve 

28/{d-e) 

^= , / V M V (11) 
(r /r) - 1 

Here r Q is the maximal radius reached by the hypersurface. 

For the case of large boundary separation L in the x direction, we can relate L to the 

depth r Q that the hypersurface reaches 

r rr r e/(d-8) 

L = drx~ / dr - = ~ r d J [d ~ e) . (12) 

J J ^{r /rf d - 1 

In particular, because 9 < d — 1 |19j . L will diverge as r a — > oo. This shows us that picking 
out the leading IR contribution, i.e. the one in the range r ~ r , is self-consistent. The 



region < r <C r Q of the integral ( 12 ) will give a contribution to L that does not diverge 
as r Q —7- oo. This is good because the IR contribution is universal in that it only depends 
on the scaling geometry ^ , while the remaining contribution is sensitive to the whole bulk 
spacetime. The case 6 — > — oo is subtle and needs to be considered separately [21| . 

We may now extract the universal IR contribution to the area and flux through the 



hypersurface. From (10), the area is 



Vol(E) 1 rd J^hT^i 



„ r d(i-d+6)/(d-e) ^ L e-(d-i) _ ( 13 ) 



While this term necessarily vanishes at large L, again because 6 < d — 1, it represents 
the leading interesting large L dependence of the area of the hypersurface. The leading 
contribution from the remainder of the spacetime will just be a constant that is independent 
of L. This constant will include the UV divergent boundary law term. The picture that 
one has in mind is that, at large separations L, the minimal surface falls straight into the 
IR region before it has the chance to significantly bend and close off. As we have discussed 



above, the behavior (13) of the hypersurface area is only sensitive to the spatial geometry 
^ and is not sensitive to the origin of the electric flux. 

To obtain the flux through the bulk minimal surface we need to know the distribution 
of bulk charges. In all known cohesive deconfined cases, where the charge extends all the 
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way into the IR, the charge density is constant in the scaling regime where the geometry 
has the form g [511 E21 122 E21 ESI ED] . This of course befits a scalar quantity in a scaling 
solution. In known partially fractionalized solutions, the charge ends a certain radius and 
does not extend all the way into the IR |29j . We consider these two cases separately. 

In the presence of a constant charge density a in the IR geometry Gauss's law for the 
(generalized if necessary) electric field E(r) takes the form 

d „ \l 9rr . \ 

—E=^—r-a. (14) 
dr r d 

Therefore in a cohesive deconfined phase, with a neutral horizon and spatial geometry Q, 
the electric field is 

oo 



E = I ^fadr ~ ar V*+0-*)/W . (15) 

J r ' 

The IR contribution to the flux through the minimal surface is then 

* r = ( Edx = r ar [d(i- d +e)+e]/(d-0) „ L 0+e/ d -(d-i) m 

Vol(S) J J r {d*-ed-B)/{d-B) ar o ■ K L0 ) 

We see that, unlike the area, the IR contribution to the flux can grow with L. Specifically, 
it grows if d(d — + 1) < 0. However, given that 6 < d— 1, it is always less than a 

volume law. 

In partially fractionalized cases the charge density does not reach the far IR. From our 



observation in equation ( 12 ), that L grows with r , it follows that the leading IR contribution 
to the flux will be insensitive to the presence of bulk charges. This is because the surface 
has 'fallen through' the region where the charges lie. Therefore Eir is constant and the flux 
through the minimal surface is 

Vq1 ^ = J Edx ~ E m L ~ (p - p coh .)L . (17) 

For the final step, we used the fact that Gauss's law implies that p = E\jy = Eir + 
Pcoh., where p C oh. is the total boundary field theory charge density that is accounted for as 
conventional gauge-invariant 'cohesive' charge. This is the charge density that contributes to 
the Luttinger count, in the case of fermionic matter. It follows that partially fractionalized 
phases exhibit the same volume law scaling for the flux <J?r as fractionalized phases, but 
with a coefficient that is less than the total charge density. 

The results we have just discussed are summarized in table [T] below. We see that the 
electric flux through minimal surfaces in the bulk is able to distinguish the four phases of 
interest. Fractionalized and partially fractionalized phases exhibit a volume law for the flux 
term, with a coefficient that detects the extent of fractionalization. Deconfined but cohesive 
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Holographic phase 


IR contribution to area Ay 


IR contribution to flux $r 


Fully fractionalized 
Partially fractionalized 
Deconfined, cohesive 
Confined, cohesive 


L e-(d-i) . Vol or L • Vol (E) 
ip-{d-x) . Yo\ (E) 
£ e-(d-i) . vol (E) 
Vol (E) 


pL-Vol(E) 

(P ~ Pcoh.) i • Vol(E) 

L e+e/d-(d-i) . Vol (s) 




Table 1: The IR contribution to the area and electric flux through a minimal spatial hy- 
persurface T in the bulk, in different holographic phases. The boundary of the hypersurface 
dT = E is two spatial hyperplanes separated by a large distance L. 



phases have a flux given by a universal L dependence that is always less than a volume law. 
Confined and cohesive phases have no flux. 

In addition to the IR contributions shown in table [T] there will be a contribution from the 
remainder of the surface. An exception to this statement is for the confined cohesive case, 
in which $r = identically. The nonuniversal contribution will be proportional to Vol (E) 
but independent of L to leading order at large separations. This provides a boundary law 
contribution. Thus we can also summarize the results as follows: 

Fractionalization =4> <I>r ~ volume law , (18) 
Deconfined cohesive =4> boundary law < < volume law, (19) 
Confined cohesive =>• <l?r = . (20) 

We emphasize that the volume law for the flux here is independent of 9 and z in the near 
horizon geometry, and in particular is not related to the volume law for the entanglement 
entropy of extremal Reissner-Nordstom black holes. 

2.2 Minimizing the deformed entropy 

Taking the deformed holographic entanglement entropy more seriously, we can select the 
bulk hypersurface V by minimizing the deformed entropy itself. Evaluating the deformed 
entropy ^ on the scaling geometry background ([9]), and using the result ( 15 ) for the electric 



field in terms of the background charge density cr, we obtain the functional 
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Minimizing the deformed entropy thus reduces to studying the dynamics given by the La- 
grangian 

C = 1 ( -L_ Vr 2e /( d - e ) + x 2 + 7 « r 9 /( d - e )^) . (22) 
r a \4Gjv / 

While we will continue to use the language of entropy, and the associated 1/(4Gjv) prefactor 
of the area, it is presumably most natural to view this Lagrangian as describing the world- 
volume dynamics of a Euclidean brane that is electrically neutral but carries an inherent 
electric dipole moment. The flux through the surface <£p ~ Jp *F precisely describes the 
coupling of a dipole moment on the surface to the external electric field. 



The quantity C = dC/dx is constant because the Lagrangian (22) does not depend on 



x(r). Knowing this, it is easy to determine that the profile of the hypersurface satisfies 

D 2 (r) 



± 2 = r 26/{d-e) g2 ( r ) ; D (r) = Cr d - 7 a r 9 '^ . (23) 



The maximal radius r a reached by the hypersurface is given by x 2 (r a ) = oo or, equivalently, 

{AG N D(r )) 2 = 1. (24) 

For 7<r = 0, this condition reduces to C = ±1/4Gtv^oi an d we recover the equation of motion 
found in the previous subsection. Thus, in particular, for fully fractionalized geometries, the 
computation and results are the same as before. The sign of C is irrelevant in this case. 
However, for 70" 7^ 0, the parameter space needs to be carefully studied. Not all values of 
7<7 need yield a single hypersurface that starts at the boundary r = and ends at r = r Q ; 
depending on the sign of C, it is possible to have none or several solutions. We will explore 



such interesting dynamical features of the Lagrangian (22) in the following subsection. 

In this subsection we wish to determine the IR contribution to the area ^4p and flux <J>p 
associated to the hypersurface T that minimizes the deformed entropy. We focus on the 
IR contribution, in the cases where this exists, by restricting all relevant integrals to the 
region r G (l — e) < r < r , with the depth r Q taken to be large. At small e, we may expand 
D(r) = D(r a ) - D'(r )(r - r), with D(r Q ) = l/AG N and 

1 



D'(r ) ~ ^ 



< 0, 

(25) 



r 



The different cases arise due to different terms balancing at leading order in the equation 
that determines the endpoint. To leading order in e, the boundary separation associated to 
depth r Q is 



L~ / drx~r e / d - e \ 

Jr (l-e) 



lr D(r 
D'(r n ) 



r d J {d - e \ 8<0. 
r (2d-6)/2( d -e) e>0 ^ 



(26) 
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In both cases L divere es as t — y oo, justifying the isolation of the near boundary contribu- 
tion to capture the large separation behavior. At this point, we are only considering values 
of Gat and 70" such that a solution to D{r ) = 1/AGm exists at large r Q for some choice of 
constant C. 

We can now evaluate the two terms in the action on these solutions. The flux is 
$r ~ / -rxr e ^ d - e > ~ { (27) 

Jr (l~e) r r {2d+e)/2(d-6)-d ^ L (0+2d(l-d+O))/(2d-9) ^ Q > q ^ 

and the area is 



r 



Jr (l-e) r 



d/(d-ff)-d ^ jl-d+9 e<Q ^ 



r o/2(d-o)-d+i _ L (-e+2d(i-d+e))/(2d-e) > fl > Q _ 

(28) 

We see that, with 9 < d—1, the area term dominates the on shell action for < 0, while the 
flux term dominates for 9 > 0. Both terms contribute equally when 9 = 0. Furthermore, 
both terms typically go to zero with large L. There is a window of positive values of 9 
just below 9 = d — 1 where the action becomes large with large L. The power of L in 
these cases is always less than one. Thus, as we concluded in the previous section for the 



flux through minimal surfaces (eq. (19)), allowing for the contribution of a non-universal 
boundary scaling term we have 

Deconfined cohesive => boundary law < Sj, < volume law . (29) 

In fractionalized phases, the absence of bulk charge combined with Gauss's law implies 
that the flux term in the deformed entropy is a total derivative that does not affect the 
dynamics. The surface T will therefore simply minimize the area. The same conclusion 
holds for the far IR of partially fractionalized backgrounds, where there is again no bulk 
charge. Our discussion in previous sections thus goes through unchanged and we conclude, 



as in eq. (18), that the deformed entropy has a volume law scaling in these cases: 

Fractionalization SZ, ~ volume law . (30) 

Finally, we consider the deformed entropy in confined cohesive phases. Let us consider 
the specific model of a fluid of charged fermions in an AdS soliton-like geometry, in which 
an internal S 1 caps off in the IR [39J. The spatial bulk geometry in the far IR is then 

^ 2 = ^ + (l-^ 2 + ^. (3D 
13 



We have rescaled the radial coordinate so that the 9 circle collapses at r = 1. The deformed 
entropy is now given by 

^ = 2vrVol(£) f dr (^L^y/ar* + U " r)x 2 + 7^=?*) . (32) 

The constant of motion is 

and the Euler- Lagrange equation can be compactly written as 



r 



2 



i „. , i ^/r^7 i 



X a i-rif 2 (r)-i ; KW_ 4Gjv r Cr d ~ l - -ytTy/Et ' ^ 

As r approaches the IR cutoff at unity, the denominator K 2 {r) — 1 becomes negative unless 
the constant C takes on the precise value C = ycry/a. Only the hypersurface with this 
constant of motion probes the deep IR; all others terminate before reaching the confinement 
scale. This unique IR-probing hypersurface has K 2 (r) diverging as 1/(1 — r) in the far IR, 
and its profile in this region satisfies 

x = ±4(d- l)a-faG N . (35) 

If this hypersurface is obtained by minimizing the undeformed entanglement entropy, 
we set 7 = and find that x = 0, namely the hypersurface falls straight into the IR 
and terminates there. This is in keeping with our previous comments concerning minimal 
surfaces in confining geometries. It is interesting that the deformation of the entropy makes 
the hypersurface profile become linear with r near the IR cutoff. This means that the entire 
bulk hypersurface gets "skewed" to one side in order to extremize the amount of flux it can 



get. Depending on the sign of ja, it is clear that one sign in (35) will have lower deformed 



entropy (32). In both cases, the minimizing surface at large separation has two disconnected 
components that fall into the far IR. Because the surfaces bend and pick up some flux in 
this case, we obtain a boundary law contribution to the flux as well as to the area of the 
surface. 

Confined cohesive ~ boundary law . (36) 

This is different to the vanishing flux we found through minimal hypersurfaces in confining 
geometries (eq. (|20[)), The deformed entropy can still distinguish confined cohesive and 
deconfined cohesive phases, however. In the confined cohesive case, there is no dependence 
at all on L at large separations, while the boundary law in deconfined cohesive phases will 
have subleading L dependence. 
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2.3 Dynamics of dipole surfaces in charged media 

Before moving on to discuss possible field theory interpretations of the above results, we will 
explore in more detail the dynamics of hypersurfaces T that minimize the deformed entropy. 
This will be the dynamics of surfaces with a tension and a dipole moment in scale-invariant 
charged media. 

We noted above that the surface has an endpoint at any radius r Q at which x(r) diverges. 



From (24), this occurs whenever 



\D{r )\ = -—; D(r)^Cr d - 7 ar e ^ d - e K (37) 



From (23) we can easily convince ourselves that this formula captures all endpoints. As we 
have already mentioned, the structure of the above equation is such that multiple endpoints 
are possible. The hypersurfaces T can consist of multiple disconnected hypersurfaces or 
of one hypersurface that does not reach all the way to the boundary. Both of these cases 
correspond to the formation of "bubbles" in the interior of the bulk. We will now describe 
the possible scenarios, with the plots in figure [4] below summarizing the various cases. 

The behavior of |-D(r)| near the boundary is primarily determined by the sign of the 
exponent 9. The fact that /{d — 6) < d for any 9 < d—1, implies that we will always have 
the scaling relation | JD (t*) | ~ r s /( d -°) at r — > 0. On the other hand, the large r behavior will 
always be |-D( r )l ~ r d . A solution x{r) will exist at those r for which |D(r )| < l/4Gjv, 
and hence we can already conclude that any solution that exists at all will have at least one 
endpoint at a sufficiently large position r Q ; there are no solutions that extend to r = oo. We 
can also conclude that, if 9 < 0, any existing hypersurfaces will have an additional endpoint 
near the boundary. Thus, no solutions reach the boundary in geometries with 9 < 0. 

Consider first 9 < 0. The function |-D(r)| decreases as we move away from the boundary, 
and at some point it will reach a minimum and then start growing again. See figure [4] below. 



Thus, eq. (37) will either never be satisfied or will be satisfied by two points r Q . In the 
former case, there is no solution. In the latter case, the hypersurface exists between the two 
endpoints, both with r a > 0, and so it will form a "bubble" in the bulk. 

If 9 = 0, the function |-D(r)| tends to \^cr\ as we approach the boundary. The minimizing 
hypersurface will reach the boundary if |7<r| < 1/4Gat. If the hypersurface does not reach 
the boundary, one may again find a bubble in the bulk. This occurs if C70" > 0; in this 
case I Z>(r) | decreases as we move away from the boundary, vanishes at Cr d = 70", and then 



increases with r until it assumes the expected IR scaling as r . Eq. (37) will then have two 
solutions corresponding to the two endpoints of the bubble. Alternatively, if C70" < (while 
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Figure 4: Typical plots of the function in various parts of parameter space. The 

horizontal lines correspond to different possible structures of solutions to |-D(r)| = 1/AGn- 



The shaded regions are where a solution to the Euler-Lagrange equation (23) exists. Bubbles 
are solutions that exist without reaching the boundary. 

Finally ifO < 6 < d—1, there will always exist a hypersurface that satisfies the equation 
of motion and reaches the boundary. If C70" < 0, the function will be monotonic, and 

it will intersect the horizontal line D = 1/4Gn exactly once, at the point r Q corresponding 
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to the single endpoint of the hypersurface. If Cjcr > 0, similar reasoning to the 9 = case 
leads to the conclusion that there are three endpoints. The solution therefore consists of 
two disconnected parts: one that reaches the boundary, and one forming a bubble. Again, 
see figure [4] 

The existence of bubble solutions may have interesting consequences for the quantum 
dynamics of our 'dipole surfaces'. It is of interest to return to these solutions once the 
microscopies of these surfaces is better understood. In terms of our immediate interest in 
solutions that reach the boundary, we see that no solutions exist at 9 < or at 9 = when 
the flux term dominates over the gravitational one, i.e. when \ ja\ > 1/AGn- In all other 
cases one can find values for the constant C for which there exist admissible solutions. In 
particular, C can be found such that the endpoint r Q becomes large. As we noted in the 
previous subsection, these solutions have the property that L = 2 xdr becomes large as 
r becomes large. Generally, there may be more than one value of C, and hence more than 
one solution, for a given large L. In this cases, the surface with lowest deformed entropy 
should be chosen 

For the cases in which there are no solutions that reach the boundary, one must consider 
more carefully the crossover between the IR scaling metrics we have been considering and 
the full spacetime, e.g. an asymptotically AdS space. In general, though, one will still find 
situations in which, at large separations L, the asymptotic solutions are repelled from the 
near horizon region, to the extent that connected solutions do not exist at large L. This 
is similar to the behavior of minimal surfaces in confining geometries. In those case one 
finds disconnected solutions at large separation. In the case of the deformed entropy in 
the presence of a background charge density, however, it is not clear what the disconnected 
solutions would be, as the surface cannot make it all the way into the IR in order to cap off 
smoothly. 

3 Towards field theoretical order parameters 

Let us discuss separately the field theory interpretation of flux through minimal surfaces 

and of the action obtained by minimizing the deformed entropy. 

Knowing the electric flux through a minimal Ryu-Takayanagi spatial hypersurface T is 

3 A detailed examination of the candidate hypersurfaces reveals that some of them may have profiles 
whose derivatives x(r) change sign before diverging. The width of such hypersurfaces can increase as we 
go further into the IR, and sometimes they can exhibit self-intersections. Our numerical studies on specific 
bulk geometries have never resulted in one such surface being the global minimizer of the deformed entropy. 
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equivalent, by Gauss's law, to knowing the electric charge contained within the surface and 
knowing the flux in the region at infinity delimited by the boundary dT = £ of the surface. 
See our discussion around figure [3] above. The flux at infinity is just the total charge of 
the dual field theory contained within the region bounded by S, which is certainly a known 
quantity. It remains to describe the bulk charge contained within the bulk surface T from a 
boundary perspective. 

The physical significance of the bulk spatial region bounded by the bulk minimal sur- 
face r is not well established. Nonetheless, it seems reasonable to expect that classical 
field configurations inside this bulk region are captured by the large N reduced density 
matrix associated to the corresponding boundary region. This is the density matrix whose 
entanglement entropy is given by the area of the bulk minimal surface. Classical field con- 
figurations should include all data that can be computed by direct evaluation on the saddle 
point 'master field' of the large N gauge theory]^] The bulk charge inside the surface is one 
such classical observable. Therefore, we may expect that the bulk charge we are interested 
in can be extracted from the large N reduced density matrix. This would show, at least in 
principle, that the flux <3?r is a well-defined field theory observable. However, it remains to 
be seen if a simple expression can be obtained for <3?r in terms of the reduced density matrix 
of the region inside E. 

We turn now to the interpretation of the deformed entropy as an action that should be 
extremized. As we noted in the previous section, this action describes a worldvolume theory 
with a tension and a dipole moment under the external electric field. To obtain the surface 
operator dual to such bulk worldvolume theories, they must be engineered (if possible) using 
D-branes in string theory. While polarization is commonplace among D-branes |43j . some 
work in a specific string theory background will be necessary to see if these effects can lead 
to the desired effective worldvolume actions. The essential question is how to maintain the 
charges on the (overall neutral) worldvolume separated, so that there is a net moment. An 
external electric field will of course polarize the worldvolume, but we are instead looking for 
an inherent dipole moment that can then interact with an external electric field. 

From an effective field theory point of view we can ask what kind of worldvolume dy- 
namics is necessary to generate or prevent the generation of a $r ~ fr*F term in the 
action. This term breaks charge conjugation C. Furthermore, the dipole moment induces a 

4 This perspective outlined in this paragraph has been emphasized to us by Matt Headrick. Related 
discussions concerning the region of the bulk captured by the reduced density matrix of a boundary region 
can be found in [ill 132] . 
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preferred orientation on the brane and therefore also breaks parity P. It follows that we can 
expect a dipole coupling to be present in the effective worldvolume action if the microscopic 
dynamics breaks C and P but preserves CP. This suggests that the surface operators we 
require could be dual to D-branes with chiral worldvolume theories. 

A third way to get a handle on the field theory physics dual to the fractionalization- 
cohesion transitions we have discussed is by generalizing the large N gauge theory approach 
to deconfinement introduced by A gauge theory on a spatial sphere undergoes a Hage- 
dorn transition at high temperatures due to the large growth in the number of long single 
trace operators. This transition describes deconfinement according to the Polyakov loop or- 
der parameter. Aspects of this story were generalized to finite density systems in [45j. The 
onset of fractionalization in large N gauge theories at finite density should be describable in 
a similar language. Instead of the partition function, one should compute the charge den- 
sity. In fractionalized phases the expectation value of the charge density will be dominated 
by long charged operators, while in a cohesive phase, the charge density will be carried 
predominantly by short charged operators. 
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